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The achievement of Bose-Einstein condensation of a dilute vapour of Rubidium 
atoms in 1995 []lj] heralded the beginning of a rapidly growing field of experimen- 
tal and theoretical endeavour, each of which has further stimulated the other. In the 
intervening years, the experimental production of Bose-Einstein condensates (BECs), 
particularly of Rubidium and Sodium vapours, has become almost routine, and atten- 
tion has shifted to either the condensation of more exotic species (of particluar note is 
the recent successful condensation of metastable helium or to deeper experimental 
investigations of the properties of an atomic BEC, up to the point of viewing a BEC 
more as a tool rather than an end in itself. The purpose here is to show links between 
chaotic dynamics, both in a classical and a quantum mechanical sense, and the possibly 
"chaotic" dynamical and stability properties of BECs. 

To begin: considering only the dynamics of a single classical point particle of mass 
m, one can reasonably speak of a class of "chaos-inducing" potentials V^(x), such that 
a particle can exhibit chaotic dynamics under the influence of such a potential. One 
can qualitatively visualize chaotic dynamics as being highly irregular motion in phase 
space. For a particle moving in one spatial dimension this is a two-dimensional space, 
where (subject to canonical changes of variables) the axes correspond to the particle's 
position x and momentum p. Such highly irregular motion is generally associated with 
exponential sensitivity to initial conditions, i.e. small differences in a point particle's 
initial position and momentum can make a large difference in its long-time behaviour. 
The appropriate classical (Hamilton's) equations of motion can then be derived from a 
Hamiltonian function of the general form: 

H = |- + K c (x). (1) 

One can readily make the jump to quantum mechanics by replacing the number- valued 
position and momentum variables with operator quantities. For a single particle, the 
dynamics are then generally most conveniently described by the following Schrodinger 
equation: 
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This is clearly a linear partial differential equation, where the time-derivative of the 
wavefunction is proportional to the Hamiltonian operator [the momenta p = 
{p x ,Py,Pz) nave been replaced by differential operators] acting on ip(x). As chaos 
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is generally associated with nonlinear differential equations, this can cause some to 
question how one can speak of quantum chaos at all. In classical mechanics we are 
interested in the motion of the particle in phase space, and thus in the changes of its 
position and momentum with time. The classical dynamics of a point particle's posi- 
tion and momentum, and the dynamics of the quantum wavefunction ^(x) are thus not 
really directly comparable: better analogies exist between Heisenberg's equations and 
Hamilton's equations of motion, or alternatively between Schrodinger's equation and 
Liouville's equation. Here quantum chaos is considered to be the study of quantum 
dynamical systems, the classical limits of which are capable of exhibiting chaotic dy- 
namics. That Schrodinger's equation is a linear partial differential equation is in this 
context an irrelevancy; nevertheless it is interesting to consider the effect of introducing 
an explicit nonlinearity into a quantum-chaotic Schrodinger equation: 
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One could regard such an equation as describing "nonlinear quantum chaos" or "non- 
linear wave chaos." Such nonlinear Schrodinger equations are relevant for nonlinear 
optics, but more pertinently are extensively used to describe BECs, when they are 
called Gross-Pitaevskii equations A Gross-Pitaevskii equation is used as an ap- 
proximate description of the dynamics of a large number N of (bosonic) interacting 
quantum mechanical particles, which is more fully described by a quantum field; if 
nearly all of these particles can be considered to be in the same motional state ?/>(x) 
(i.e. if we have a BEC), then it is possible to describe most of the underlying quan- 
tum field by a classical field, and to regard what is left over as being "small." The 
Gross-Pitaevskii equation is then the equation of motion of this classical field. It is 
nevertheless sometimes necessary to consider the dynamics of the full quantum field; 
in second-quantized form, the Hamiltonian operator leading to such a Gross-Pitaevskii 
equation (with u = gN and g being a parameter essentially describing how strongly 
individual particles interact) is given by 



H = J d 3 xv;t (x) + y c(x ) + |^t (x) ^ (x) 



$(x), (4) 



where V^( x ) an ^ VH X ) are bosonic field operators, which create and annihilate, re- 
spectively, a particle at position x. In this context, one could thus possibly speak of 
"quantum field chaos" when considering the dynamics induced by such a Hamiltonian. 

A brief discussion of integrability, a concept closely tied to the (non) chaoticity of 
dynamical systems, is now necessary [Q]. Recalling the classical dynamics of a sin- 
gle particle, such a system is considered integrable if there exist as many independent 
conserved quantities as there are motional degrees of freedom (i.e. the number of rel- 
evant spatial dimensions). If the system is integrable, then the phase-space dynamics 
of the particle are so restricted by the necessity of observing these conservation laws 
that chaos is impossible. To take the rather trivial example of the one-dimensional har- 
monic oscillator [potential V(x) — mui 2 x 2 /2], the solution of Hamilton's equations 
of motion (in slightly unconventional form) is given by: 

a(t) = e -<w *a(0), (5) 
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where a is a complex variable, defined as 

mux + ip 

a= -=-. (6) 

V2m 

It is clear that the particle motion is restricted to closed circles in the two-dimensional 
phase space, making both "irregular" motion in phase space, and sensitivity to initial 
conditions completely impossible. In the language of Hamilton-Jacobi theory, \a(t)\ 
is an action variable, which is conserved (and equal to the square root of the energy), 
and the dynamics are fully described by the time-evolution of the phase angle of a(t), 
which is then the canonically conjugate angle variable 6(t) to |a(t)|, defined by 



V mux / 



arctan 
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Similar analyses can be made for any one-dimensional, time independent system; due 
to the time-independence of the Hamiltonian, the energy is always a conserved quan- 
tity, which is enough to make a one-dimensional system integrable. Two- and three- 
dimensional systems will require additional independent conserved quantities to be in- 
tegrable, and time dependence in a one-dimensional system generally lifts the property 
of integrability. Concepts of quantum integrability generally, by analogy, consider con- 
served observable quantities, which in quantum mechanics are described by operators 
[Q]; it is possible to consider something different however. 

A wavefunction tp(x) can generally be decomposed in terms of some complete 
orthonormal basis [St]: 

oo 

VK X ) = ^2c k ip k {x). (8) 

k=Q 

If additionally the ip k (x) are eigenfunctions of the (assumed time-independent) Hamil- 
tonian, then any time dependence of ip(x) is contained completely within the coeffi- 
cients c k : thus 

c fe (i) = e-^*c fe (0). (9) 

For example, in the case of a one-dimensional quantum harmonic oscillator, the ip k (x) 
are Gauss-Hermite polynomials, and u k = u(k + 1/2). There is an now an obvi- 
ous similarity between equations (|j) and (Q). Extending the analogy, we can consider 
motion in an infinite-dimensional pseudo phase space of a "point particle" having "po- 
sition" values 

Xk = c k +c* k , (10) 

and "momentum" values of 

Pk = ~i(c k ~ c%). (11) 

Within this pseudo phase space however, as in the case of the classical harmonic oscil- 
lator, motion is severely restricted by the fact that all the \c k \ are conserved The 
\c k \ (effectively the "action variables" in the pseudo phase space) in fact form a com- 
plete set of conserved quantities (i.e. there are as many conserved quantities as degrees 
of freedom), meaning that, in this sense, the Schrodinger equation is always integrable 
and thus never chaotic, a fact which is indeed a direct consequence of the Schrodinger 
equation's linearity. 
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For a general Gross-Pitaevskii equation, the additional nonlinear term means that 
integrability clearly cannot be taken for granted. It is, for example, known that the one 
dimensional, homogeneous case 
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is integrable am, and that the conserved quantities can be determined using inverse scat- 
tering techniques. However, even for the seemingly trivial extension of simply adding 
a harmonic potential (generally more relevant to current experiments), it is unknown 
whether the resulting Gross-Pitaevskii equation is integrable or not; the application of 
a 5-kicking potential (as is widely considered in the study of classical and quantum 
chaos) 



V c (x) = K cos(kx) 



S(t-nT), 



(13) 



will produce a Gross-Pitaevskii equation which is unlikely to be integrable [g], This 
general lack of integrability therefore raises the possibility of chaotic dynamics in the 
pseudo phase space, including extreme sensitivity to initial conditions. 

In classical machanics, sensitivity to initial conditions is generally gauged by the 
Lyapunov exponent, so that a positive Lyapunov is considered a strong indicator of 
chaotic dynamics. The Lyapunov exponent is concerned with the growth with time of 
the euclidean "distance" d in phase space between two initially close point particles. 
In one dimension, d is thus given by 



d(t) = y/[ Xl (t)-X 2 (t)]^ + \ Pl (t)- P2 (tW 



(14) 



In the pseudo phase space produced by some orthonormal representation of two initial 
wavefunctions ipi(x) and tp2{x), the equivalent d is given by 



2- / d 3 x[^(x,^2(x,t)+V 2 *(x,i)V'i(x,i)], (15) 



where in the case of the linear Schrodinger equation, by unitarity J d 3 x-0i (x)^2(x) 
cannot change, and d will therefore be constant |^]. Upon the addition of a nonlin- 
earity this is of course no longer guaranteed, which in the case of the Gross-Pitaevskii 
equation can have important consequences for the stability of the condensate, due to 
a correspondance between the propagation of linearized perturbations around the con- 
densate wavefuntion, and propagation (and growth) of the non-condensate fraction. 

It is now necessary to return briefly to a fuller quantum mechanical description 
of the BEC. If the number of particles in the condensate mode is nearly N, the field 
operator f/>(x) can be written as [|Io[ |IJ 



V-(x) 



V'(x) 



kfcWfc (x) + b\v* k (x) 



(16) 
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The a operator annihilates a particle in the condensate mode ^(x) (it is approximately 
this mode which is propagated by the Gross-Pitaevskii equation if the number of par- 
ticles in the condensate mode a) a w AO, and the St and are creation and annihi- 
lation operators of excitations above the condensate, the (x) and vt (x) describing 
the spatial dependence of these excitations. These excitation modes clearly describe 
non-condensate particles, and for the idealized case of zero temperature, the number 
of non-condensate particles is then given by 

J d 3 xk(x,i)| 2 . (17) 

k 

The excitation modes are propagated by slightly modified Bogoliubov equations, which 
are identical to the equations derived by considering perturbations around and orthog- 
onal to ^>(x) propagated by the Gross-Pitaevskii equation up to linear order in the 
perturbation only. If one then considers 

V>i(x,0)=V(x,0), (18) 

and 

^ 2 (x,0)=V(x,0) + ^=[u fc (x,0)+^(x,0)], (19) 
v N 

it is clear that a reduction with time of | J d 3 xipl (x, t)ip2 (x, t) | directly implies growth 
in the number of non-condensate particles, and, more generally, growth in the number 
of non-condensate particles is implied by any kind of linear instability in a solution 
of the Gross-Pitaevskii equation. Chaotic dynamics in the pseudo phase space thus 
imply the rapid increase in number of non-condensate particles, which in turn implies 
rapid depletion of the condensate. There is a caveat to this: the derivation of the equa- 
tions used in the analysis formally take the limit N — ► oo, g — ► 0, while keeping 
gN constant. There is therefore formally an infinite number of condensate particles, 
which stays infinite as the number of non-condensate particles increases, and the Gross- 
Pitaevskii equation is always valid. 

In any real situation the total number of particles must of course be finite. Recall- 
ing that it is assumed that nearly all the particles are in the condensate mode ^(x), 
it is clear that once depletion starts to become significant, this approach can be prob- 
lematical. In the systematic (and equivalent) expansions of C.W. Gardiner jlQ], and of 
Castin and Dum JTT[], there are in principle further terms which can be incorporated to 
account for higher order effects. Morgan JI^ ] has taken a slightly different approach 
by only requiring that the occupation of the condensate mode be large, i.e. it may be 
significantly different from the total particle number, and adding in perturbative cor- 
rections to the Gross-Pitaevskii and Bogoliubov equations. This treatment, however, is 
designed for studying essentially steady-state BECs at finite temperature, rather than 
the explicitly dynamical sorts of situations described here. 

The theoretical treatment of such situations is clearly a difficult problem, it is how- 
ever likely relevant to coming generations of BEC experiments. If one actually wants 
to do something with a BEC, this will involve dynamics. As generic Gross-Pitaevskii 
equations are presumably non-integrable, there exists the definite possibility of chaos 
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and instability in some circumstances; this in turn implies the possibility of rapid de- 
pletion of the condensate mode, and concomitant loss of the desirable property of co- 
herence of the macroscopic matter wave. Given this, it is equally desirable that such 
processes be better understood; this thus seems a suitable "quantum challenge" for the 
21st century! 
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